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SUMMARY 

A brief resume of sonic boom minimization methods is given to provide a background for 
a new, empirical modification of the Seebass and George minimum-nose-shock sonic boom F- 
function and signature. The new “hybrid” F-function has all the inherent flexibility of application 
found with the Darden-modified Seebass and George F-function. In addition, it has enhanced this 
flexibility and applicability with neglegible increase in nose and/or tail shock strength. A 
description of this “hybrid” F-function and signature is provided, and the benefits of using them to 
design high-performance, low-boom aircraft are discussed. 


INTRODUCTION 

The theory of minimum sonic boom has advanced considerably since the L. B. Jones report 
on a lower bound to sonic boom, reference 1, in 1961. Based on the flow field of a supersonic 
projectile paper by G. B. Whitham, reference 2, Jones predicted that a slender body whose area 
growth was proportional to x 1 ^ would produce a minimum disturbance overpressure signature in 
the far field. 

Walkden, reference 3, showed that the lift on a wing-body also contributed to the far-field 
overpressure because, along the longitudinal or flight direction, it generated equivalent area which 
was directly proportional to the product of the local lift and the Mach number parameter, p, and 
was inversely proportional to the cruise velocity dynamic pressure. With this extension to 
Whitham’s theory, the L. B. Jones lower bound body could now represent both the volume and the 
lift of the aircraft However, the nose bluntness of the lower bound body meant that aircraft it 
represented would be subject to a sizeable zero-lift wave drag penalty. 

While evaluating the overpressures generated by large aircraft during the acceleration to 
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cruise velocity, reference 4, F. E. McLean noted that equivalent-area body shapes with less wave 
drag than the Jones lower bound body produced lower intensity shocks at distances in the mid-field 
range. Studies of sonic boom propagation through a standard stratified atmosphere by W. D. Hayes, 
reference 5, demonstrated that the shape of the pressure signature would tend to stabilize or 
“freeze” while the shock strengths attenuated as the disturbances travelled the mid-field distances 
between the stratospheric cruise altitude and the ground. These studies changed the focus of 
minimization from the far-field to the mid-field. 

The low boom and minimum boom signatures used today came from the mid-field studies 
of R. Seebass and A. R. George, reference 6, who applied minimization techniques to pressure 
signature shapes. They provided a means by which the overall aircraft design could be guided and 
controlled to reduce the shock strengths felt on the ground. These methods still led to aircraft with 
some low-boom-induced drag penalties although they were usually lower than previous levels 
found on the lower bound equivalent area body. However, with the modification to the nose- 
bluntness requirement introduced by C. M. Darden in reference 7, trade-offs between shock 
strength and aircraft drag could be conducted during the design phases that showed promise of 
providing aerodynamically efficient, mission-capable aircraft with acceptable sonic boom 
characteristics. 

Another modification to the overpressure signature shape recently contributed by G. T. 
Haglund, reference 8, further opened the “window” of design flexibility for the aircraft designer 
integrating both low sonic boom and high aerodynamic efficiency characteristics into the aircraft 
configuration. Like the Darden nose-bluntness modification, it was suggested by a merging of 
experience with the Seebass and George minimum boom pressure signatures with purely practical 
design considerations. The purpose of this report is to show why and how this “second generation” 
sonic boom signature was developed and demonstrate the benefits derived from applying it to 
conceptual aircraft design methodology. 


SYMBOLS 

A e aircraft equivalent areas, ft 2 

B value of the F-function slope between y = ^ and y = 1 , ft 

v 

C value of the F-function between y = yf and y=^,ft 
D discontinuous change in the F-function at y = X, ft 

F(y) the Whitham F-function 

h cruise altitude, ft 

H value of the F-function at y = yf / 2 

I(x) unit step function; 1(a) equals 1.0 only for x ^ a , it is zero elsewhere 
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K f ground reflection factor, usually 1.9 

^ effective length of the aircraft or model, ft 

M cruise Mach number 
p flow field pressure, psf 

p a ambient pressure, psf 

Ap p - p a > Psf 

q free stream dynamic pressure, psf 

W beginning cruise weight, lb 

x longitudinal distance, ft 

Xp effective distance, ft 

y effective distance parameter in the F-fimction, F(y) 
yf F-function “nose-bluntness” parameter, ft 

y r effective distance along trailing F-function where tail shock conditions are met 

a value of the atmospheric “advance factor'’, AF(y)/Ay 
a atmosphe ric advance defined in reference 6, a = C / a , ft 

p Jm 2 -i.o 

X effective length that determines the shape of the positive portion the F-function, ft 

£ effective length where “ramp” of slope B begins, ft 

<t> angle whose tangent is equal to the value of B, see figure 7 


DEVELOPMENT 

L. B. Jones, one of the first to research minimum sonic boom body shapes, derived his 
lower bound body from impulse theory. It has rapid local area growth at the nose followed by 
smaller area growth such that the equivalent area distribution of this lower bound sonic boom body 
is proportional to x 1/2 . With this F-function, an N-wave shaped far-field overpressure signature 
with the lowest ground overpressure level is obtained. However, the drag increment associated 
with the local nose bluntness can be prohibitively large. Figure 1 shows the equivalent area 
distribution, the F-function, and the pressure signature of a typical Jones lower bound body. 
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Figure 1 . Area distribution, F-function, and pressure signature of a Jones lower bound body. 


When later research and design studies indicated that supersonic-cruise aircraft would be 
long and slender enough so that the mid-field rather than the far-field signature would reach the 
ground, other types of body shapes were examined for both low boom and low wave drag 
characteristics. The minimization techniques used by Seebass and George provided two pressure 
signatures which were constrainted for either minimum nose shock or minimum peak 
oveipressure. Both of these are shown in figure 2 as presented in their original form. 
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Minimum Overpressure 


Minimum Nose Shock 



Figure 2. Seebass and George minimized pressure signatures, equivalent areas, and F- 
functions. 


Note that both of these signatures come from F-functions with Dirac delta-functions at their 
origins, just like the Jones lower-bound-body F-function. They are different in that the distributed 
area behind the delta-function permits a lower nose shock solution for a specified cruise altitude 
than that provided by the Jones lower bound body for an aircraft of the same length and weight. 

By replacing the Dirac delta function with a narrow adjustable nose “spike”, as shown in 
figure 3 for the minimum overpressure signature F-function, the drag of the vehicle can be reduced 
at the cost of a small increase in shock strength. 
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With this modification, the configuration drag can be varied during the design stages. 
Trade-offs between drag and ground overpressure permit more flexibility in selecting and 
integrating aircraft components. 

Although not shown, the modified Seebass and George minimum shock or “ramp” F- 
function, reference 6, can be obtained by the addition of a “spike”. In addition to the nose “spike” 
width, yf, the slope of the “ramp” length, B (see symbol list and figure 2), can also be adjusted to 
permit more component arrangement flexibility in the aircraft design. 

Both of the minimum boom F-functions and signatures are point-design shapes with point- 
design limitations. The “flat-top” signature has one forward of the y =yf point on the F-function, 
while the “ramp” F-function and pressure signature has one on each side of the y = yf point. With 
the “ramp” F-function, perturbations to the ambient conditions of the standard atmosphere will 
change the atmospheric propagation characteristics resulting in higher overpressures. The “flat- 
top” F-function is somewhat less sensitive in that atmospheric perturbations in only one direction 
will produce higher overpressures. 

These design point features indicate that the “ramp” and the “flat-top” signatures can be 
minimized only for a narrow range of atmospheric propagation characteristics. This poses no 
problem when designing research or wind tunnel models, but is a potential difficulty when 
designing a real supersonic cruise aircraft. Since the standard atmosphere is an averaged statistical 
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model, aircraft configurations designed with it and these two F-functions have limited value 
because only an average set of flight conditions are being considered and met. 

G. Haglund made a further practical modification to the F-function and pressure signature 
to overcome these point-design limitations. Like the nose “spike” feature which preceded it, the 
modification was empirical and increased the applicability and flexibility of the low-boom 
methodology at the potential expense of a small increase in nose shock strength. The derivation of 
the equivalent areas from the F-function and a description of the F-function features are outlined 
in the following section. 

In figure 4, this “second generation’-modified F-function, pressure signature, and 
equivalent area distribution are presented. 



Figure 4. Haglund-modified equivalent areas, F-function, and pressure signature. 

This F-funcdon is no longer a true minimized F-function nor does it give a true minimized 
pressure signature in the classical mathematical definition. However, for practical aerodynamic and 
acoustic reasons, it combines the best features of the Seebass and George F-functions and the nose 
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“spike”(for drag- nose shock trades), with a constant-value section between the “spike” and the 
“ramp”. Such an F-function and its derivative signature is more accurately called a “hybrid”. 

The constant-value section aft of yf in the F-function permits hot-day and cold-day 
perturbations in the atmospheric propagation characteristics to be incorporated into the aircraft 
design as well as allowing additional control over aircraft length, component arrangement and 
integration, and area growth. With this feature on the F-function and its incorporation in the aircraft 
geometry, the probability that the ground shock strength would be predictable for a specified Mach 
number, altitude and beginning cruise weight could be greatly improved. The following paragraphs 
will describe these two features. 

Atmospheric compensation. The method for incorporating hot-day, cold-day atmospheric 
propagation variations can be explained with the use of figure 5. 


Advance lines 



Figure 5. Hybrid F-function showing hot-day/ cold-day compensation analysis. 


For given flight conditions of cruise Mach number, altitude, ground overpressure, and 
weight, the atmospheric propagation advance is computed for the standard day, the hot day, and the 
cold day. Cold-day conditions decrease the value of the atmospheric “advance factor” resulting in 
a spike length that is longer than one required for a standard day. Conversely, hot-day conditions 
increase the atmospheric “advance factor” value relative to that for a standard day and the “spike” 
length is shorter. This inverse relationship between the atmospheric “advance factor” and “spike” 
length is due to the definition of the “advance factor” which represents the net change in a finite- 
pressure signal location relative to a zero-strength acoustic signal traveling through the 
atmosphere. 

The lower value of the advance factor could determine the location of the “nose bluntness” 
length, yj , while the highest value of the advance could determine the most forward location of 
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jjje “ramp” starting point, £ . With these two values in hand, the remaining signature parameters of 
Mach number, cruise altitude, aircraft effective length, slope of “ramp”, and cruise weight would 
^ | n put to the Hybrid Signature code. Computed values of Ap would be compared with desired 
overpressure limits. Adjustments to the input, excluding y f and £, would be tried iteratively until 
satisfactory F-function, equivalent area distribution, and signature was found with % no less that 
the value set by cold-day conditions. 

Usually, the “spike” length, yf , is varied to study the trade-off between shock strength and 
zero-lift wave drag. It still can be used for this purpose as long as the hot-day/cold-day atmospheric 
advance compensations are not compromised. There can still be length available for a larger value 
of ^ which would be useful in giving sufficient aircraft volume to meet mission requirements of 
fuel volume, passenger cabin room, reserve fuel, landing gear stowage, etc. Two examples are 
given to show the sensitivity of the method. 

The first is for an aircraft which cruises at a Mach number of 3.0 : 

M = 3.0 , 
h = 73,000 ft, 

W = 650,000 lb, 
l e = 300.0 ft , 

% = 40.0 ft , and 
B = 0.5 a 

For a 1962 standard atmosphere day, the value yf = 20.0 ft results in a Ap = 1 .032, psf. 
To obtain the same value of Ap on a “hot day”, yf needs to be about 18.0 ft, while for a “cold day”, 
yf would be about 22.0 ft. So if an overpressure of 1.032 psf is an acceptable nose shock strength, 
then 18.0 feet value provides the necessary “cold day” compensation. The aircraft will have a bit 
more nose bluntness and probably a bit more drag than if it were designed for a standard day, but 
it will meet the desired nose shock criteria for all but the most extreme atmospheric conditions. The 
hot-day and the cold-day atmospheres were approximations to those defined in the 1962 standard 
atmosphere tables. However, the results obtained with their use indicated the probable values and 
the ranges in “hybrid” signature-calculation constants determined with more exact methods. 

The second example is for an aircraft which cruises at a Mach number of 1.6 : 

M= 1.6, 

h = 45,000 ft , 

W = 650,000 lb , 

l e = 300.0 ft , 

£ = 40.0 ft , and 

B = 0.5 a 
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For the standard day, y f = 20.0 ft results in Ap = 0.617 psf. For the “hot day”, yf of about 
19.0 ft is appropriate, while for the “cold day”, y f is about 2 1 .0 ft. The plus-and-minus differential 
of “hot day” minus “cold day” y f values are very small at either Mach number. Thus, only the “hot 
day” value is of primary importance. Most of the time, the value of £ is greater than the yf value 
required for “cold day” compensation. This permits additional flexibility in adjusting the desired 
low boom equivalent areas for good agreement with those of the conceptual aircraft being 
designed. 

The Ap value of 0.617 psf looks attractive but should not be taken seriously without 
further examination. In this example, the aircraft nose geometry which would produce such a 
ground nose shock strength was prohibitively slender. Thus, nose geometry as well as overpressure 
is a consideration in the selection of an equivalent area distribution. 

Design flexibility. The advantages of using the Haglund “hybrid” F-function method during 
the preliminary design phases are illustrated in figure 6. 



In this example, the Mach 1.6 configuration, previously mentioned, serves to demonstrate 
the increased design flexibility acruing from the use of the “hybrid” F-function derived equivalent 
area distributions. Its beginning cruise specifications are repeated for ease of referral. 
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M = 1-6 . 

h = 45,000 ft , 

W s=650,000 lb , 

i =300.0 ft , 

c 

£ = 40.0 ft , and 

B = 0.5 a 

A value of about 19.0 feet for y f will meet “hot day” requirements but results in an aircraft 
nose that is unusually slender. Decreasing yf to 10.0 feet while keeping the rest of the input values 
constant makes the nose blunter and structurally more practical. It also reduces the overpressure. 
However, this nose shock overpressure reduction will come with the addition of more aircraft 
volume, and potentially more drag, even though the weight is assumed to have remained constant. 

Two pairs of lines are shown in figure 6. The upper pair is for the total equivalent area due 
to lift and volume, while the lower pair is for the fuselage areas. Reducing yf increases the aircraft 
volume. This volume increment can be added to the fuselage to permit six-abreast rather than five* 
abreast seating. If this is not necessary, the extra volume can be used to obtain additional wing 
volume for fuel. Often, drag and/or weight penalties result from increasing the volume while 
keeping the length constant. If the drag and weight increments from this volume increase are found 
to be relatively small, this new low-boom solution may be more suitable than the previous one. The 
equivalent area c'urves shown in figure 6 for values of y f equal to 10.0 and 19.0 feet indicate that 
over a range of about 220.0 feet in effective length, an increment of about 10.0 square feet in 
equivalent area would be required. 

With the exception of the fuselage, this example is based mostly on overall equivalent areas 
rather than on actual geometry. If an actual conceptual aircraft were being designed, it is very 
possible that some of the other parameters would also have to be changed to obtain a good 
agreement between the aircraft and the ideal equivalent area curves. Obviously, the shape and 
location of the aircraft components are also changing, but by varying both the aircraft geometry 
and the low boom equivalent curve, the time required for convergence can be reduced. This double- 
effort approach assumes that a specified overpressure level is never compromised. 

Assuming that the Mach number, the beginning cruise altitude, and the beginning cruise 
weight are fixed, the values of B and ^ remain as variables to be altered as the design matures. In 
contrast, the modified Seebass and George F-functions permit changes only in the value of B after 
fixing the size of yf and maintaining the other parameters as constants. In the next section, the 
derivation of the equivalent areas equation is presented and described. The conditions for 
determining the strength and position of the nose and trailing shocks are also described to show 
how they are used to determine the necessary constants and coefficients in the area equation. 
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DERIVATION OF THE “HYBRID” EQUIVALENT AREAS 

The F-function shown in figure 4 is inverted by using Abel’s integral in the form 
x 

A e (x) = 4.0 JF (y) Jx-y dy 
0 

using the same constraints as in reference 6. It is repeated as figure 7 so the various features 
can be explained. 



The equivalent area associated with this “hybrid” F-function is: 

Ae(x) = (32/15) [(H/y f ) x 5/2 + I(y f /2) [x-jj +1 (y f ) 

(x-y f ) 5/2 + I(^)(l/2)B(x-^) 5/2 -I (X) (5/4)D(x-X) 3/2 ] 
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where I(x) is the unit step function. 


At -x = l e 

A e (le) = (P/2)W/q 

plus any equivalent area adjustments for wake, boundary layer displacement thickness, 
engine-nacelle exit area minus intake area increments, etc. Nose shock Ap is computed from the 
F-function parameter, C , the “advance factor”, a , the altitude, h , the reflection factor, K r , and 
the ratio of cruise altitude and ground pressures. 


As it is applied in this paper 


a = AF(y)/Ay 


This “advance factor” should not be confused with the “advance” used in the reference 6. 
There, the advance is called a x , the distance that a unit disturbance would lead or follow an 
acoustic signal that reached the ground after propagating through the atmosphere. It would be equal 
to 


a x = C/a 


The “advance factor” is the change in distance with altitude that a unit strength disturbance 
travels relative to an acoustic signal starting from the same point along the aircraft during cruise 
flight. Like the advance, the “advance factor” is calculated from the Mach number, the altitude, and 
the characteristics of a standard, “hot day”, or “cold day” atmosphere. The “advance factor”, a , 
relates the two F-function parameters H and C through 


H = 


C 2 

y f a 


C 

2 



C 

2 


The triangular “spike” of magnitude H at y = yf/2 and C at y = yf , is the modification 
introduced by Darden in reference 7 to reduce the nose bluntness associated with the areas derived 
from the Dirac-delta function on the F-functions described in reference 6. At y = £ , the constant 
value of F(y) = C , the Haglund innovation, ends, and F(y) continues with slope B past the 
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discontinuity at y = £ to y = 1 . The value of AF(y) = D at y = X .the slope B , and the 
aircraft or body length 4 are used to set the ratio of tail shock strength to nose shock strength. A 
solution for the tail shock is found from the value of F(1 ), a , and the integral of F(y) between 1 

V V 

and an F-function area-balancing point yr which is solved iteratively through 

'r 

0.5 (F (y r ) - F (l e ) ) (y r -l e ) = J F (y) dy 
and 


a(y f -l e ) =F(y r )-F(l e ) 

Values of yf , l e , B , and £, are part of the input parameter set. Other input parameters 
include Mach number, altitude, ground reflection factor, and ratio of tail shock strength to nose 
shock strength. The nose shock strength, Ap , is an output value computed from the input values 
and shock conditions; it is not a specified or target input. If it is not satisfactory, some of the input 
lengths, the altitude, the weight, or the Mach number will have to be changed. 

The code that computes the hybrid F-function, equivalent areas, and signature is meant to 
supplement the Seebass and George minimization code. Together, they allow the designer to obtain 
the possible minimums and then trade sonic boom and atmospheric perturbations with aircraft drag 
and performance to obtain a satisfactory sonic boom constrained configuration. . 


CONCLUSIONS 

Past experience with sonic boom minimization methods and techniques have shown that 
the pure mathematical approach has produced two types of performance penalties. The first type 
was associated with zero-lift wave drag due to the locally severe blunting applied to the nose of 
the aircraft. This was found with the L.B. Jones lower bound sonic boom body and to a lesser 
extent, the Seebass and George set of minimum nose shock and minimum overpressure signatures. 
The second type of penalty was due to the point-design nature of the F-functions or the 
corresponding pressure signatures. Predicted overpressure signatures might be obtained for a 
specified standard, hot-day, or cold-day atmosphere, but all three conditions could not be satisfied 
simultaneously. 

A set of empirical cures were found to overcome these limitations. Using a triangular 
“spike” rather than a Dirac delta-function permitted drag-shock strength trade-offs to be studied 
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and applied. Combining the “flat-top” and the “ramp” F-function shapes and starting the 
p-function with the previously-mentioned nose “spike” provided the additional flexibility 
necessary to meet drag constraints as well as variable atmosphere ambient conditions. The 
resulting “hybrid” F-function and pressure signature was not a minimum in the mathematical 
sense, but was a practical compromise in terms of the airplane configuration itself. 

Hot day-cold day conditions are small perturbations to the standard day features of the 
Hybrid F-function. While providing a useful limit on the nose-bluntness length, yf , they should 
not seriously hamper efforts at setting the “ramp” initial length, £ , and the “ramp” slope, B, such 
that aircraft volume, aircraft mission performance, and low sonic boom constraints can be met. 
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